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ABSTRACT 


With the use of a diffuse potential of the Woods-Saxon type, including a spin-orbit interaction, 
numerical solution of the single nucleon Schrédinger equation is performed, with values of the 
potential parameters appropriate to the lead region. Single particle energies and radial wave 
' functions are obtained for the orbitals in the major shells adjacent to Z = 82 and N = 126, and 
the changes of the energy eigenvalues with the parameters are estimated. Pairing energy correc- 
tions are introduced. The possible occurrence of quadrupole oscillation states is considered, and 
the corresponding collective corrections to the single particle states are applied. The perturbed 
energies are compared with the level schemes of the four nuclei, which differ from the double- 
magic Pb®°8 with one nucleon, and some suggestions and predictions are made. 


1. Introduction 


The use of a nuclear potential with a diffuse surface in shell model calculations 
has been necessitated by experience [1]. In the optical model analysis of cross sections 
for the elastic scattering of nucleons from nuclei, a potential of the form 


V(r) = —V, [1 + exp(r — R)/a} 


introduced by Woods and Saxon [2], has often been used [3, 4]. Calculations of bound- 
state energies with such a potential plus a spin-orbit interaction with the radial 


dependence 
ridV/dr 


were first performed by Ross et al. [5]. Recently, Sliv and Volchok [6] have studied 
the nuclear potential parameters with the same assumptions. Similar calculations, 
though with a somewhat different potential form, have been performed by, among 
others, Bleuler and collaborators [7]. 

One of the main objects of the earlier investigations has been to obtain a proper 
choice of the potential parameters common to the whole periodic system. A more 
detailed analysis for some specific regions, where the shell model is expected to 
hold, might be of interest. In the present work, the single particle states in the region 
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2. The nuclear potential and the eigenstates “a F 
2.1. The potential and the eigenvalue problem 


- According to the model applied, the interaction of a single nucleon with the residual 
nucleus can be represented by three contributions to the Hamiltonian: A spin-inde- 
pendent central nuclear potential, a spin-orbit coupling term, and, in the case of a 
proton, a repulsive electrostatic potential, here referred to as the Coulomb-term. 
We assume the following forms: 

a) The central nuclear potential 


V(t) Survnmameety regi te elie ARh: (1) 


V, > 0 is the depth of the potential and R its radius. With R =7,A"®, ry is approxi- — 
mately constant over the whole periodic system. 2a is a measure of the depth of the 
surface. 

b) The spin-orbit term 
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A is a positive dimensionless parameter, characterizing the strength of the coupling. — 
A, = h/Mc is the reduced Compton wave-length of the nucleon. The angular momentum _ 
operators § and / are given in units of hi. 

ce) The Coulomb term 
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(Z—1)—E(r), (3). 


where Z is the proton number and £(r) is chosen to correspond to a uniform charge 
distribution within the radius R 


? 


546 


a Hein, MW? 
avast % ela hs “ 
orien ee ter 


ote 
Paging 


, sy y iny tt nae ee ei an Qu ian Deerpety 
pepe OHA Forage Oe amnge fra MY RUGS, 
; 1) for j=l—1/2, age 


n is the usual [1] radial quantum number (n- 1 equal to the number of 
Ray (r) for 0<r< oo), 
‘ insertion of the function 


Yny (7) =7 Ray (7) (6) 


ae 


and of numerical values for the constants (units: MeV and 10°¥ cm), the 
equation reads 

_ ah = ” L(i+1 
a ¥()- 


y (r) + 0.04823 F + 


+Vo E (r) + 0.02209 f (, j) an) y (1) — 


( 
. 
q 
: 
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— 0.06939 (Z—1) ie GE (7) 


_ The subscripts nlj have been omitted for simplicity. The boundary and nor- 
_ malization conditions are 


y (0) =y (ce) =0, a ET ae, a) 


We put Z — 1 =82 in the following. 

The unphysical singularity of &’(r)/r for r = 0 is removed by setting &(r)=1 ina 
small region r < b< R. The appropriate solution y(r) is then proportional to r times 
a spherical Bessel function for r< b; see equation (A.4) in Appendix A. For the 
normalization one needs a fairly correct representation of y(r) for large values of 
r,r —R>a. The asymptotic form is given in Appendix A 

The eigenvalue problem (7), (8) has been treated numerically with the Swedish 
electronic computer BESK. Given the potential parameters, the eigenvalue EH for 
a state n1j can easily be estimated with sufficient accuracy for the adopted proce- 
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dure. Equation (7) is then integrated with the method of Runge-Kutta, starting at — 


the point r = (a suitable choice of b is 0.5- 10-1 cm). The following rule holds: 


(—1)"y(r)> +00 if B>Ey, 


t> oo 


(-1)"y(r)>-—c if B< Ey, 


to 


where E,y is the exact eigenvalue. The difference H — H,,, varies smoothly with 
y(r) at a sufficiently large value of r. These facts furnish a rapidly convergent itera- 
tive procedure for getting the eigenvalues very accurately. The asymptotic form of 
the eigensolution y(r), however, has to be calculated by other means, and the 
expression (A.6) has been used for this purpose. The “tail” is then fitted to the 
integrated solution y(r) at a suitable point r=d (d has here generally been chosen 
as 13-10-13 cm). Finally, the functions y(r) have been normalized according to (8). 


2.2. Parameters and results 


The choice of the nuclear potential parameters is given in the first row of Table 1. 
They are compared with the values obtained by other authors from shell model caleu- 
lations [5, 6] and optical model analysis [8, 4]. 

The number of decimal figures, which we have chosen for our parameters, is such 
that another decimal would certainly be physically insignificant. The parameters 
have not been fitted in detail to any specific level data, but have been chosen to fulfil 
the following approximate conditions, which fix their values fairly well. 

(I) The values should be in rough agreement with the results obtained by other 
authors for the whole periodic system, according to Table 1. We set the limits: 
R = (7.6 + 0.3)-10-18 om, V2 =45+5 MeV, V2? =57+5 MeV, a = (0.7 +0.1)- 10-8 
om, 4=35 +5. 


Table 1. The potential parameters. 


R Vo MeV 
Reference : OG A aie "i A 
10-13 ¢ 10-13 @ —13 
wa Rae Neutron | Proton alin 
The present choice for the 
Pb-region 1.273 7.62 44.0 58.0 0.67 32 
Pipe? gt 42.6 — 0.87? 28 
Sliv and Ppt? é 7.691 41.8 — Orr 39 
; 4 cone 1.3 i ; 
Volchok [6] Bi? %. 71+ ~— 56.9 0.87? 25 
AWEW 7.691 — 55.8 0.73? 38 
Ross et al. [5] 1.3 er Oe 42.8 56.6 0.69? 39.5 
Elastic scattering of 
neutrons [3] 1.25 7.411 50° — 0.65 358 
Elastic scattering of 
protons [4] 1.25 7.41} — 60° 0.65 358 
: ones appropriate to the Pb-region according to the radius law R = 1 Alls, 
“a«=a" given in the references. 


° Rough estimates, 
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to the potential parameters. 


Table 3. The partial derivatives of the energy eigenvalues with respect 


; Peeinced, as it containe the best eee date in the region (00 Appendix 0). 
_ Especially the difference between the energies of the states 1)2 (ground state) and 
—%132 (experimentally 1.63 MeV) gives a sensitive criterion for the choice of A. 


\ 


Neutron states 


State oH oH MeV aH MeV OF 2 Mev 
nj aVo éR 10-8 cm da 10-8 em i 

W desis —0.44 -4.9 aot + 0.022 
2 Guya ~0.59 Oy al ~0.8 + 0.052 
4 84/5 — 0.42 =43 =32 0 

Ted ~0.77 =8.9 +7.6 —0.107 
Sides — 0.55 —5.5 ~2.5 ~ 0.016 
Weis -0.73 ~9.2 +707 + 0.086 
2 Gojo —0.69 — 6.8 +07 — 0.041 
3 pas — 0.69 -6.4 -0.1 + 0.019 
2 fsje — 0.73 ==i7ieb) +2.2 + 0.041 
3 Paya ~0.72 ~6.2 +0.3 —0.009 
Visa ~0.81 —7,3 +8.0 — 0.084 
2 faa — 0.78 ~6.4 +2.8 —0.031 
IV — 0.80 18 + 8.5 + 0.063 

wee la ee ee ee ee EE eee 

Weta — 0.83 —10.0 +113 — 0.109 
2 frye OSL ors + 4.5 — 0.039 
1 fgje — 0.82 --11.0 siz ek’) + 0.079 
eas —0.83 67 AL 0 

2 dyin — 0.84 Hot +5.8 + 0.034 
lh A — 0.86 = fey) + 10.9 — 0.080 
2 ders ~ 0.86 —8.4 +5.7 — 0.023 
eee = 0°87 -9.4 +115 + 0.056 


Proton states 
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Fig. la and 1b. The normalized function Yn1y (7) 
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i normalized eigenfunctions y(r) are tabulated in Appendix B and plotted in 
Figures la-1d. The accuracy of the calculations is such that we expect the error to 
be at most one unit in the last decimal place. The relatively long “tails” in some of 
the cases of Fig. 1a depend on the low energy values in these cases. 

The dependence of the energy eigenvalues on the nuclear parameters is shown in 
Table 3. The partial derivatives 0#/0R, 0E/0Vo, 6H/éa, and 6H/éA have been 
evaluated numerically, using formulas of first order perturbation theory. In some 
cases the results have been tested by integration of (7), (8) for slightly different values 
of the parameters. The derivatives are in fair agreement with values obtained by Sliv 
and Volchok [6], the deviations being due to the different parameter values used. 

‘Since shell model calculations-have often been made with the use of harmonic 
oscillator radial wave functions (cf., for instance, True and Ford [8]), it is of some 
interest to expand the more physical wave functions, which have been obtained with 
a diffuse potential, in these oscillator functions. This has not been done in the present 
work. Sliv and Volchok [6] have investigated this point, and their results should 
hold essentially unchanged for the functions obtained here. 


3. Corrections and comparison with data 


As is well-known, there are—due to residual interactions—different correlations 
between the shell model particles. These turn up more or less pronounced in all parts 
of the periodic system, and are essentially of two types: The low order correlations, 
which give rise to pairing effects, and those of high orders, which come into evidence 
in collective vibrations and rotations. Even for the low-lying states of the four nuclei 
adjacent to Pb’, for which we are going to make a detailed comparison of theoretical 
and experimental energies, there are effects of the correlations, which cannot be 
neglected. For the lowest levels, the energy corrections are not too large, and this 
makes it possible to caleulate them approximately. Some simplifying assumptions 
are necessary in order to avoid involved calculations and the introduction of many 
new parameters. The change of the corrections due to variations of the potential 
parameters can be neglected, which makes the procedure of choosing parameters 
much easier. Special attention is paid to the possible occurrence of quadrupole vibra- 
tional states, which cannot be simply described in terms of single particle excitations. 

The consistency of the method depends on how adequately the different effects 
are represented (single particle states; pairing energies; collective vibrations). One 
may hope to get a better fit to the data with the corrections applied than with the 
pure single particle states, 
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a. The pairing energy corrections 


The interaction causing pairing effects is here represented by a 6-function. The 


pairing energy for an orbital n/j is then to the first order given by the diagonal 
- matrix element 


6= —K(294+1)L,y= 


foe} 


a 
3 = —K(27+1){ [Ray (nr)? dr. (9) 
%, 0 


_ The quantities 6 are added as corrections to the single particle energies corresponding 
to states in the nuclei Pb?°’ minus one nucleon. We neglect the possible effects of 
_ scattering of pairs between the major shells. The value K = 6-10-89 MeV- cm? is 
chosen for both the neutron and proton cases, giving reasonable agreement with 
experimentally estimated pairing energies. The overlap integrals J,,,, and the correc- 
tions 6 are given in Table 4. 
As the admixtures to the wave functions due to pairing effects are negligible for 
our purposes, we will continue to speak of single particle states, even after the pairing 
corrections have been applied. 


b. The collective corrections 


An odd nucleon or nucleon hole, moving in the field of Pb?°, has a polarizing in- 

. fluence on the otherwise spherically symmetric, closed shell nucleus. This is a non- 
static effect which quantum-mechanically can be described as a collective oscillation 
of the nuclear surface coupled to the single particle motion by an interaction operator 
Hj. We consider only harmonic quadrupole vibrations. In the general coupling 
scheme, the states are then denoted |/j; N R; J M>, where N is the number of vibra- 
tion quanta, phonons, and F their resulting angular momentum. Furthermore, we 
will use the one-phonon-approximation, i.e. disregard the effects of states with 


Table 4. Pairing energy corrections. 


Overlap ‘pn ete 
State integral box: 3 poe 
nlZ Ur 
1 030" ern MeV 
Neutron hole states 3 Pre 0.0137 —0.16 
(corresponding to Pb?%’) 2 fsje 0.0084 — 0.30 
3 P3/2 0.0141 — 0.34 
1 43/2 0.0068 — 0.57 
2 fajs 0.0083 — 0.40 
1 hg/s 0.0082 — 0.49 
Proton hole states 3 81/2 0.0392 — 0.47 
(corresponding to Tl’) 2 ds/p 0.0124 — 0.30 
1 Ayo 0.0082 — 0.59 
2 dsia 0.0119 — (0.43 
1 97/2 0.0091 — 0.44 
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Table 5a. Collective corrections to the single-particle states in the cases 
Pb208 + one nucleon. Some of the lowest one-phonon states. 
hw =2.5 MeV. Admixture: Squared amplitude. 


i TE | 


; : . A rrected 
Predominant Main admixture Correction re 
Nucleus pe configuration oe BiNhote 
l ] NR °, i i N’ R’ e Mev 
Ns ee SS SS eee 
Pb?” 9/2+ aan 00 2. eens 12 -0.11 — 3.80 
11/2* Hag! 00 2 Aa 12 ~0.11 3.07 
5/2+ ee 22 Joja 12 —0.42 — 2.24 
15/2- Liziteat 00 3 iis 12 —0.07 i 
1/2+ aun m0 8 ae 12 —0.27 —1.49 
7/2+ gr ~~ 55 90/2 12 —0.43 eA, 
3/2+ 3/2 00 4 I7/2 12 — 0.28 — 0.88 
7/2* Ioj2 12 37 Ii 00 —0.11 ~1.3 
11/2+ Io/2 12 0 0.00 —1.2 
9/2+ EYP) 12 2 Jo/2 00 + 0.06 —1.1 
5/2+ eaten tal? 20 pee 00 + 0.29 ~0.9 
Bi209 9/2- Fan 00 2 Thess 12 — 0.07 — 3.61 
t[2- 7/2 00 Be frie 12 — 0.07 — 2.76 
13/2+ ead 2 ile 12 —0.07 — 2.34 
Ue hoje ie 0 + 0.01 —1.0 
9/2- i), eae 2 hea 00 +0.07 LtKO 


1 This configuration is not predominant, but is obtained in the limit when the coupling para- 
meter eé is 0. 


more than one phonon. According to Bohr and Mottelson [9], the relevant matrix 
elements are given by! 


<19; 00; 7M | Hine |l’9’; 125 7 MY =F e<j|h|7">, (10) 


where ¢? = e- hw, with iw = the phonon energy and e = k?/2C (k is the coupling con- 
stant and C the deformability). ¢j|h|j’> are real numbers, for which formulas are given 
in reference [9]. In all four cases we put e = 0.69 MeV, in agreement with an estimate 
by True and Ford [8] from the probability of the 0.57 MeV £2 transition in Pb207. 
This value is probably of the right order of magnitude for nuclei adjacent to Pb, 
k is essentially constant, about 40 MeV [9]. C on the other hand, varies strongly 
and becomes large in the double-magic regions. The quantitative information is very 
incomplete in the lead region, but we assume C to have a value between 1000 and 1500 
MeV, as has been found for Pb2°? and Pb [8, 10]. 

The evidence for fiw is also weak. The energy of the first excited 2+ state for 
spherical even-even nuclei is known to increase rapidly as one moves towards closed 
shells (reference [10], especially Fig. V.3). In Pb’, no 2+ state has been found, but 
there is no strong evidence against it above 2.5 MeV [8]. One might expect that ha 
is of the order of 2 MeV for the nuclei Pb2°8 + one nucleon. We will use Aw = 1.5 MeV 


1 — for particle, + for hole. 
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+ 0.16 


+0.21 
+ 0.25 
+ 0.06 
—0.07 
—0.00 
—0.09 
=u.1T 
+ 0.08 
+ 0.03 
| — 0.06 
T1207 1/2+ uae 00 4 Can 12 +0.15 — 7.35 
3/2+ ae 00 4 ‘Sine 12 +0.16 — 7.65 
11/2- ky 00 4 hap le + 0.06 ~9.19 
5/2+ dex 0 9 ae 12 —0.07 — 9.84 
7/2+ Iza 00 15 ee 12 — 0.03 STs) 
5/2+ eS 12 dea 00 Sen —8.9 
3/2+ ee tb 4 djs 12 — 0.05 =9.1 
7/2+ Soa 22 2 Guia 00 + 0.05 = 9 9 
1/2+ tg wie 4 ore 00 —0.08 aos 


i 


for Pb2°? and T1297, and fiw = 2.5 MeV for Pb? and Bi?. The former choice is made 
on the following grounds. With reasonable values of the potential parameters, the 
fr2 neutron hole level in Pb?’ always comes out 0.5-1.0 MeV higher in energy than 
the experimental 7/2- level at 2.35 MeV (cf. reference [6]). We suggest that the latter 
is mainly an f5)2 + one phonon state. This assumption, which will be discussed later, 
then fixes the value of fo. 

The calculation of the collective energy corrections A is performed by a straightfor- 
ward diagonalization of the energy matrix, where the non-diagonal matrix elements 
are given by (10) (if not zero), and the diagonal ones are E + Nio for the particle 
cases, and H +6 —Nfw for the hole cases. This procedure also gives the perturbed 
wave functions. The results are given in Tables 5a and 5b, where we have also 
included some of the lowest perturbed one-phonon states. 

A survey of experimental data for the nuclei Pb?, Bi, Pb0?, and T1?°7 is given 
in Appendix C. The theoretical and experimental energies are compared in Figures 
2a and 2b for neutron and proton states respectively. Those levels, which have been 
suggested to be of core-excited (c.e.) or octupole vibrational (0.v.) character, have 
no counterparts among the theoretical levels in Figs. 2a and 2b. 

As far as level data are well-known, the agreement between the calculated and the 
experimental energies is satisfactory. This fact supports the consistency of the method. 
It may therefore be justified to go somewhat further in the interpretation of experi- 


mental data and to make some predictions. 
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Fig. 2a. Theoretical unpertubed (to the left) and perturbed energies for neutron states, and 
experimental levels of the nuclei Pb?°* + one neutron. Dotted lines denote the lowest quadrupole 
vibrational states. Arrows indicate energies used for the choice of parameters. ‘‘c.e.”” means 
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_ Fig. 2b. Theoretical unperturbed (to the left) and perturbed energies for proton states, and 
a experimental levels of the nuclei Pb?° + one proton. Dotted lines denote the lowest quadrupole 
_ vibrational states. Arrows indicate energies used for the choice of parameters. “‘o.v.”” means 
a : octupole vibrational state. 
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the vibrational g,/. + one phonon state (22% with ou 


meters). . -_ ae 
ioe 2.13 MeV level is probably a core-excited (4/2) (Yo/2)” state aS 
suggested by Stephens [12]. The energy agrees well with the difference between the 
neutron binding energies for Pb®°8 and Pb”!®, 2.14 MeV. - 

Our results point to an s,/. rather than a ds). assignment for the 2.01 MeV level. 
The energy of the theoretical s,;. state could be lowered by using a potential form 
with a somewhat slower decrease for r >-R. : 

Pb7,—In choosing potential parameters, the positions of the p,/. and tg/2 States 
have been adjusted to agree roughly with those of the corresponding experimental ‘ 
levels. The fs. and pg. states are seen to fit fairly well into the scheme. 

The experimental 7/2- level at 2.35 MeV should have mainly the structure of an 
fs neutron hole coupled to one phonon. There are two arguments which support 
this interpretation. In the first place, the 1.78 MeV y-transition between the 7/2~ and 
5/2- states goes by ~ 94% H2 and ~ 6% M1as deduced by Alburger and Sunyar [13] 
from a y —y angular correlation measurement.? If this were a single neutron transition, 
one would expect the multipolarity M1 with a very small admixture of #2. Secondly, 
the 1.78 MeV transition is stronger than the 1.44 MeV transition to the 3/2- level by 
a factor of 50, and hence the ratio between the reduced £2 transition probabilities 
is ~16. Both these facts indicate a collective 7/2-—5/2- transition. 

The levels at 2.71 MeV and higher energies, which have been seen in the Pb? 


(d,p) Pb? reaction [11, 14], are probably core-excited states with a structure like © 


the ground state of Pb°* with one neutron added in the next major shell. The f,,. and 
hg neutron hole states are not expected to show up in this reaction due to the small 
admixtures of (f,/.)~? and (hg/.)~? in the Pb°* ground state. 


{ 
' 
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The possible vibrational 5/2- and 3/2- states occurring at an excitation energy of — 


~ 1.7 MeV would not be populated with sufficiently large intensities to be seen either 
in the (d, p) reaction or in the 6 decay of Bi?®’. 

The 9/2- state from the coupling of an f,,. neutron hole to one phonon would be 
reached by electron capture from Bi?°? with an intensity comparable to that of the 
7/2- state from the same configuration if it occurs at the same energy or lower. We 
have actually found the 9/2~ state to come somewhat lower than the 7/2- state, but 
perturbations from states with two phonons may shift their order, as illustrated in 
intermediate coupling calculations by Choudhury [15]. The present experimental 
information does not exclude the existence of a 9/2- state just above the 7/2- state. 

Bi? .—For the first excited state, we have regarded the assignment 7/2- as estab- 
lished. The theoretical results predict the assignment i,,/. for the second excited state. 
The psjp, f5/2, and py. states have been found to be unbound, but a potential form 
with a slower decrease for r > R might bind the p,/. and fs). States. Virtual single 
proton states should, for instance, give rise to low energy proton groups in the 
Pb?°8 (d,np) Pb298 reaction. 

The experimental level (or levels) at 2.6 MeV, which is strongly excited in inelastic 


' Note added in proof: This is also suggested by G. Holm, J. Burwell, and D. Miller from 
results of (d, p)-reactions on Pb?°8 and Bi? (to be published in Phys. Rev.). 
* Note added in proof: Predominant M1 character has, however, been obtained from polarization - 


direction measurements by P. H. Stelson, W. G. Smith, and F. K. McGowan (to be published 
in Phys. Rev.). 
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it can be considered as a core-excited state. . 
t be that it belongs to the lowest group of quadrupole vibrational 


- 


x ’ : : 
The energy separation of the two lowest states is well reproduced by the 
2 I results, and the suggested spin and parity assignments are confirmed. The 
olute energies, however, disagree by about 0.7 MeV. The agreement can be con- é 
erably improved by changing the potential parameters for protons, while keeping 4 
1em unchanged for neutrons, but this has not been done as the significance of 
such a procedure is doubtful. “00” 
The h,,/, state is expected to decay by an M4 transition to the dj). state, unlike 
the lighter odd Tl-isotopes, for which the h,,). state should decay with multipolarity 
#3 to the experimentally observed 5/2+ state at 0.6-0.8 MeV [17]. Because of the 
_ small separation of this 5/2+ state and the d,. state, the former can hardly be inter- 
_ preted. as ds ;., but should rather be due to neutron excitation and should therefore 
not occur in T]?°7, 


4, Discussion 


: For the potential terms (1), (2), (3) and the involved parameters we have made the 
_ following basic assumptions: 
1. The nuclear potential V(r) is a monotonic function, and its variation is sym- 
metric around the point r = R. 
- 2. The parameter V, is energy independent but has unequal values for neutrons 
- and protons. 
3. Except for the difference between Vj and V¢, the potential parameters are the 
same for the different nuclei. 
4, The spin orbit term has the radial dependence r+ V’(r). 
5. The Coulomb term is chosen to correspond to a uniform charge distribution. 
The last assumption is not quite correct. One should rather use a charge density 
with a smooth decrease at the nuclear surface. The details of the charge distribution 
are, however, not crucial for the purpose of the present work. Assumption 4 is some- 
what arbitrary but physically reasonable, and the results are not very sensitive to 
the radial dependence of the spin-orbit interaction, provided it is concentrated to 
the nuclear surface. 
The condition 1 is fulfilled by the mathematically convenient form (1) of the 
diffuse potential. As mentioned above, one might in some cases get a better fit to 
_ the data by using a potential with a slight asymmetry around r= Rk. This could be 
achieved by “filling up” the function V(r) in a suitable interval tong aah 
Assumption 2, about the energy independence of Vo, is not altogether satisfactory. 
From optical model analysis [3, 4], V_ is found to decrease with increasing nucleon 
bombarding energies, and one might expect the same tendency to hold for bound 
states. This is also one of the results of Brueckner’s theory of nuclear matter [18]. 
However, it seems as yet difficult to draw any definite conclusions about the magni- 
tude of the variation of V, with energy for bound states in a nucleus. The main effect 


1 This modification would also reduce the separation of the perturbed d3/2 and h11)2 levels to 
about 1.0 MeV, which indicates that we may have overestimated the energy difference between 


the 3/2+ and 11/2- states. 
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“core” of 126 neutrons and 82 protons. Thus the radius # has been kept at a constant 
value appropriate to A = 208, and the factor Z—1 in the Coulomb term (3) has been 


put equal to 82 for both the proton and the proton hole cases. Accordingly, no at- — 


tempt has been made to trace a dependence of V4,a, and A on the individual nucleus, 
as has been done by Sliv and Volchok [6]. In this respect our approach is essentially 
different from theirs, as we apply corrections due to collective and pairing effects. 
Yet the nuclear parameters may have some individual variations, which we have 
neglected. 

The justification for the choice of the collective parameters iw and e has been 
discussed above. The values of iw chosen for A = 207 and A = 209 are rather uncer- 
tain but not physically unreasonable. e may vary somewhat from nucleus to nucleus, 
but due to the lack of information we have put it constant, equal to the value ob- 
tained for Pb?°’. 

The choice of the parameter K in the expression (9) for the pairing energy may be 


subjected to criticism. Matrix elements for the short range forces used by True and — 


Ford [8], in good agreement with experimental data for Pb®°*, correspond to a value 
of K about twice as large as our choice. This comparison is, however, not very con- 
clusive. If a larger K-value is assumed, one must also, to some extent, adjust the 
potential parameters. We have not entered upon this readjustment, particularly as 
it should not essentially influence the results or any conclusions drawn from them. 

Finally one might discuss the conditions (I)-(II1) for the choice of the potential 
parameters, especially concerning the use of the nucleon binding energies. However, 
the details of the assumptions are not very critical. The requirement (III), which may 
seem to be somewhat arbitrary, helps to fix the parameters V¢, Vj, and R and gives 
the approximately correct binding power of the nuclear potential. The separations of 
the levels within the major shells, on the other hand, furnish tests for the parameters 
a and A. With the data available so far, it is probably not possible to determine the 
potential parameters in an essentially more consistent way. 

In conclusion one would like to express the expectation that many new experiments 
will be done to investigate the properties of the low-lying levels of the nuclei Pb297, 
Pb??, T?°7 and Bi2. The experimental information is as yet generally too scanty to 
admit more detailed theoretical investigation along the present lines. 
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In the same interval, the function ¢, in (3a) can also be given a non-essential change so 


iva 


ee: calth = oR constant. 
_ The equation then becomes, for r<b, 
ee - ” L(l+1) 
: were |e°— ae] y (r)=0, (A.1) 
_ with 
j B=(0—xe-BY, (A.2) 
x = 0.04823 V5; k? = 0.04823 | EL bs 
be. (A.3) 
L Tl 
; x2 = 0.10409 a (0 for neutrons). 
- 
4 The solution with y(0)=0 is 
2 y (1) = Nye Tis, (BP), (A.4) 


4 where N, is a normalization constant. With the actual choice of b, the Bessel-function 
can be sufficiently well approximated by the first few terms of its series expansion in 
powers of #7. From this one can also easily calculate the logarithmic derivative of 
y(r) at the point r = 6, which is used in the numerical treatment of the eigenvalue 


problem. 
a The asymptotic form of the wave function 
From a point r = d and outwards, the function & (r) and its derivative are negligible, 
and the wave equation thus reads, for r> d, 
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A certain approximation method has been used for computing this function in 
intervals of interest. In the case =0 (neutrons), the Whittaker function redu 

a Hankel function with imaginary argument. For 7 >/ there are special difficult eS. 
However, it is possible in all the cases to compute the “tails” of the wave functions 
with sufficient accuracy for the present purposes. 


APPENDIX B jy 
The functions y(r) according to (6), 


y(r) =r R(r), 


where K(r) is the normalized radial wave function, are given in Tables B 1-B3. 


N 


APPENDIX C 
Experimental data 


Of the four “‘single-particle nuclei”, Pb2°7 is stable, Bi2® stable or nearly so, and 
Pb? and TI’ are 6--emitters. From the tables of Wapstra [19], one finds the neu- 
tron binding energies of Pb®°® and Pb28 and the proton binding energies of Bi? 
and Pb’ to be 3.87, 7.38, 3.72, and 8.05 MeV respectively. Experimental excitation 
energies, spins, parities, and decay properties of levels in the four nuclei are given 
in Table C 1. Assignments, which are considered tentative or not quite certain, have 
been put in parenthesis. 

The experimental results for Pb2 are rather scanty. Some information has been 
obtained from (d,p)-reactions [11, 14], some from spectroscopical investigation [12]. 


The assignments are not known with certainty, not even for the ground state, although 
9/2+ seems very probable. 
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7 453 -508 | 4420) "672 -504 | 453 


7.5 477 — 457 477 | 469 -476 | 334_ 
8 428 — 369 . | 432 347 -391 | 230 
8.5 950 1" Wa. = 278 353 240 —293 152 
9 272 — 202. 272 159 — 209 98 
9.5 206 145 203 103 — 146 63s 
10 154 — 108 150 67 —~ 100 a OT 
10.5 115 =" 109 44 —68 26 | 
11 85 - =52 80 29 —47 1st 4 
11.5 63 — 37 58 19 — 32 ll 
12 47 ~26 42 12 ~22 7 
12.5 35 -19 31 8 -15 5 
13 26 =48 23 6 $0 3 
13.5 19 -—10 16 4 -7 2 
14 14 =9 11 3 -5 1 
14.5 Ua) -—5 8 2 -3 1 
15 8 =4 6 1 = 1 
15.5 6 -3 4 1 -2 0 
16 4 -2 3 1 -l 
16.5 33 -] 2 0 -l1 
17 2 -1 2 ai 
17.5 iP. -l 1 0 
18 1 -]l 1 
18.5 1 0 l 
19 1 1 
19.5 1 0 
20 0 


The positions of some low-lying states of Bi2 are well established by various 
experiments, mainly inelastic nucleon scattering [16, 20]. The ground state spin is 
known from atomic spectroscopy to be 9/2. 

In Pb? four excited states below 2.5 MeV are well-known from spectroscopical 
investigations of the decay of Bi297, especially the work by Alburger and Sunyar [13]. 
Harvey [11] obtains from (d, p)- and (d,t)-reactions all the states given in Table C 1, 
and also further states up to 5 MeV. More recent verifications are furnished by the 
results from the (d,p)-reactions of McEllistrem et al. [14] and the (p,p’)-reactions 
of Cohen and Rubin [16]. Atomic spectroscopy gives the spin 1/2 of the ground state. 
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The information about Tl2°? is very meagre, and only one y-ray has been attributed 
to that nucleus [21], supplementing the information from «-spectroscopy. Systematics 
of odd TI! isotopes indicate the ground state spin 1 /2. Cf. Bergstrém and Andersson 
uae 
Concerning the ground state moments of the four nuclei, the situation is briefly 
the following. For Pb?°® and TI? nothing is known experimentally. The magnetic 
moment of Bi2® is +4.1 nuclear magnetons, and the deviation from the Schmidt 
value +2.6 has been attributed to configuration mixing. For Pb?” the experimental 
value agrees with the Schmidt value, +0.6 nuclear magnetons. The quadrupole mo- 
ment of Bi2°° is —0.3 (in units of e-10-*4 cm?), and the shell model value is —0.2. 
A possible difference may be ascribed to weak surface coupling, but gives no sig- 


nificant information about the collective parameters. 
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1.634 13/2+ | 0.8 sec, ie 64. Mev, M4 
es ieee MeV, E2(~94%) +. 

2.35 ! 1.45 MeV, y J 

2.71 (9/2+) a ie 

3.61 | @A/2+) 

T1207 0 (1/2+) B- 
0.350 (3/2+) 0.350 MeV, M1+E2 
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